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Within an isospin- and momentum-dependent hadronic transport model it is shown that the 
recent FOPI data on the n~ /n + ratio in central heavy-ion collisions at SIS/GSI energies (Willy 
Reisdorl et al, NPA 781, 459 (2007)) provide circumstantial evidence suggesting a rather soft 
nuclear symmetry energy E sym (p) at p > 2po compared to the Akmal-Pandharipande-Ravenhall 
prediction. Some astrophysical implications and the need for further experimental confirmations 
are discussed. 
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The density dependence of the nuclear symmetry en- 
ergy E sym {p) is critical for understanding not only the 
structure of rare isotopes [l| and heavy- ion reactions @, 3, 
3, la , but also many interesting issues in astrophysics 6, 



To determine the E sym {p) and thus the equa- 



tion of state (EOS) of neutron-rich nuclear matter has 
been a longstanding goal of both nuclear physics and 
astrophysics. While both fields have some promising 
tools for probing the E sym {p) over a broad density range, 
they all have some limitations. Thus, only by combin- 
ing carefully complementary information from both fields 
will we have ultimately a good understanding about the 
E sym (p) ■ While significant progress has been made over 
the last few years in constraining the E sym (p) at sub- 
saturation densities using terrestrial nuclear laboratory 
data@, still very little is known about the E sym (p) at 
supra-saturation densities. 

In fact, the high density behavior of the E sym (p) has 
long been regarded as the most uncertain property of 
dense neutron-rich nuclear matter [13, EH- Presently, 
at supra-saturation densities even the trend of the 
E sym {p) , i-c, whether it increases continuously or de- 
creases at some point with the increasing density is 
still controversial. While many microscopic and/or phe- 
nomenological many-body theories using various inter- 
actions, such as the Relativistic Mean Field [l2j and 
Brueckner-Hartree-Fock approaches [IH, predict that 
the E sym (p) increases continuously at all densities, an 
approximately equal number of other models includ- 
ing the Variational Many-Body theory [3, EH, EH] and 
the Dirac-Brueckner-Hartree-FockE3 predict that the 



and the density-dependent M3Y interactions [24|, 125J or 
the non-relativistic Thomas-Fermi approach using the 
Seyler-Blanchard interaction [26| also lead to decreas- 
ing E sym (p) starting at some supra-saturation densi- 
ties. Thus, currently the theoretical predictions on the 
E S ym{p) at supra-saturation densities are extremely di- 
verse. Therefore, to make further progress in deter- 
mining the E sym (p) at supra-saturation densities, what 
is most critically needed is some guidance from exper- 
iments. In this Letter, we report circumstantial ev- 
idence suggesting a rather soft nuclear symmetry en- 
ergy E sym {p) at p > 2po compared to the Akmal- 
Pandharipande-Ravenhall (APR) prediction [53] based on 
a transport model (IBUU04[28]) analysis of the recent 
FOPI data on the tt~ /tt + ratio from central heavy- ion 
collisions at SIS/GSI energies [H. 

The isospin and momentum-dependent mean field po- 
tential (MDI) used in the IBUU04 reads HI 
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where p n and p p denote the neutron (r = 1/2) and pro- 
ton (r = —1/2) density, respectively, and 5 — {p n — 
Pp)/(Pn + Pp) is the isospin asymmetry of the nuclear 
medium. All parameters in the above equation can 
be found in refs.pci]. The variable x is introduced to 
E sym (p) first increases to a maximum and then may mimic different forms of the E sym (p) predicted by vari- 



start decreasing at certain supra-saturation densities de- 
pending on the interactions used [HI. Additionally, the 
non-relativistic Hartree-Fock (HF) approach using many 
Skyrme [H, El, M, E3 , Gogny [H], Myers-Swiatecki [23| 
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ous many-body theories without changing any property 
of the symmetric nuclear matter and the E sym (p) at nor- 
mal density p . The E sym (p) with x values of 1,0.5,0 
and —1 arc shown in Figure CD By setting x = 1 
one recovers the HF prediction using the original Gogny 
force [30l|. As it is well known, the latter predicts a 
soft E sym(p) decreasing with increasing density at supra- 
saturation densities. For comparisons, shown also are the 
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E sym {p) = 12-5p/po + 12.7(p/p ) 2/3 used in the IQMD 
(Isospin-Dependent Quantum Molecular Dynamics) HU 
and the well-known APR prediction [27]]. The latter has 
been widely used in calibrating other model calcula- 
tions. In terms of reproducing the experimental data, 
the IBUU04 model has had modest successes so farQ. 
While the NSCL/MSU isospin diffusion data[H allowed 
us to limit the E sym (p) at sub-saturation densities to be 
between that with x = and x = — 1 [33l |34| , the same 
model parameter sets underpredict 13511 significantly the 
double neutron/proton ratio of ref.[3y]. Nevertheless, it 
is interesting to mention that the above limited range 
of the E aym (p) for p < po is consistent with that ex- 
tracted very recently from analyses using the ImQMD 
(Improved QMD) model which can reproduce both the 
isospin diffusion data and the double neutron/proton ra- 
tio simultaneously [37| . It is also worth noting that the 
APR prediction for the E sym (p) at sub-saturation densi- 
ties lies right between that with x = and x = — 1. 
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FIG. 1: (Color online) Density dependence of nuclear symme- 
try energy E 3y m(p) predicted by APR (stars), used in IQMD 
(dotted line) and that using the MDI interaction with x=l, 
0, 0.5 and -1, respectively. 

Among the most sensitive probes of the E sym (p) at 
supra-saturation densities proposed in the literature 
the 7r~/7r + ratio in heavy-ion collisions is particularly 
promising. Qualitatively, the advantage of using the 
7r~ /n + ratio is evident within both the A(1232) reso- 
nance modeljll] and the statistical model[3^| for pion 
production. Assuming only first chance inelastic nucleon- 
nucleon collisions produce pions and neglecting their re- 
absorptions, the A resonance model predicts a primordial 
7T-/7T+ ratio of (7r-/7r + ) ros = (5 N 2 + NZ)/(5Z 2 + NZ) « 
(N/Z) 2 cnsc , where the N and Z are neutron and pro- 
ton numbers in the participant region of the reaction. 
The 7r~/7r + ratio is thus a direct measure of the isospin 
asymmetry {N / Z)& cnac of the dense matter formed. The 
latter is determined by the E sym (p) through the dy- 
namical isospin fractionation [Jjj, namely, the high(low) 
density region is more neutron-rich (poor) with a lower 

Es y m 

(p) at supra-saturation densities. Since effects of 
the E S y m (p) are obtained mainly through the correspond- 
ing nuclear mean-field which dominates the dynamics of 
heavy- ion reactions at relatively low energies, based on 



the resonance model one thus expects the it /ir + probe 
to be most effective at beam energies near the pion pro- 
duction threshold E£ h w 300 McV. On the other hand, 
assuming pions have gone through multiple production- 
reabsorption cycles and reached thermal-chemical equi- 
librium, the statistical model predicts that 7r~/7r + oc 
exp [2{p n - p p )/T] = exp [8SE sym (p)/T], where T is the 
temperature. Thus, in this model the n~ /-k + ratio mea- 
sures directly the E sym (p) at the pion freeze-out. Mean- 
while, at energies much higher than the E£ h where pi- 
ons are abundant, the reaction dynamics is dominated 
by scatterings among all hadrons instead of the nuclear 
mean- field 4l| . Therefore, one expects that the tt~ /it + 
probe becomes less effective at very high energies where 
other observables, such as, the neutron-proton differen- 
tial flow[42|, EH, may be more useful for probing the 
high density E sym (p) . More quantitatively and real- 
istically compared to the above two idealized models, 
several hadronic transport models have shown consis- 
tently that the ir~ /ir + ratio is indeed sensitive to the 
E S ym (p) [IS H El especially near the E£ h . Moreover, 
by varying separately the E sym (p) at sub- and supra- 
saturation densities in IBUU04 simulations we found 
that the n~ /n + ratio in collisions near the E£ h is much 
more sensitive to the variation of the E sym (p) at supra- 
saturation rather than sub-saturation densities. 

Recently, Reisdorf et al. studied systematically the 
Ca+ 40 Ca, 96 Ru+ 96 Ru, 96 Zr+ 96 Zr 



r /tt + ratio 



40/ 



and 197 Au+ 197 Au reactions using the FOPI detector at 
SIS/GSl[29|]. Their Tr~/-K + data are the most exten- 
sive and accurate ones available in the literature, thus 
providing us the best opportunity so far to extract the 
E sym (p) at supra-saturation densities. 
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FIG. 2: (Color online)The ir~ /tt + ratio as a function of the 
neutron/proton ratio of the reaction system at 0.4 AGeV with 
the reduced impact parameter of b/b max < 0.15. The inset is 
the impact parameter dependence of the tt~ /tt + ratio for the 
96 Ru+ 96 Ru reaction at 0.4 AGeV. 

Shown in Fig. [5] are the calculated tt~ /tt + ratios in 
comparison with the FOPI data at 0.4 AGeV with the 
reduced impact parameter bo = b/b max < 0.15 as a func- 
tion of the neutron/proton ratio of the reaction system. 
The inset shows the ir~ /tt + ratio as a function of bo for 
the 96 Ru+ 96 Ru reaction at 0.4 AGeV. It is seen that 
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both the data and the calculations exhibit very weak bo 
dependence for the tt~ /?t + ratio, even for mid-central re- 
actions where we found that the multiplicities of both ir~ 
and 7r + vary appreciably with the bo . For the symmet- 
ric 40 Ca+ 40 Ca and the slightly asymmetric 96 Ru+ 96 Ru 
reactions, calculations using both x = 1 and x = can 
well reproduce the FOPI data. Interestingly, for the more 
neutron-rich reactions of 96 Zr+ 96 Zr and 197 Au+ 197 Au 
calculations with x = 0, 0.5 and 1 are clearly separated 
from each other. The FOPI data favors clearly the cal- 
culation with x = 1. As shown in Fig. [TJ with x = 1 
the E sym (p) at p > 2po reached in the reaction is very 
small, leading to a rather high N/Z in the participant re- 
gion and thus the larger 7r~ /tt + ratio observed. We note 
here that in the present study we used directly the con- 
dition bo < 0.15 in selecting the most central events. On 
the other hand, in the data analyses and the IQMD cal- 
culations, the distribution of ERAT (ratio of transverse 
to longitudinal kinetic energies) was used in determining 
the centrality of the reaction. Nevertheless, given the fact 
that the tt~ /it + ratio is almost a constant within error 
bars over a large range of bo, more elaborate selection of 
the most central events by matching the calculated and 
the experimental ERAT distributions may change some- 
what quantitatively the 7r~/7r + ratio, but it is not ex- 
pected to change qualitatively our conclusions. For com- 
parisons, the IQMD result from ref.[29[ is also shown. As 
seen in Fig. [Q the E sym (p) used in the IQMD, the MDI 
E sym (p) with x=0 and the APR prediction are all very 
close to each other for po < p < 3po- Thus, not surpris- 
ingly, the 7T-/7T+ ratios from the IQMD and the IBUU04 
with x=0 are very close too. Moreover, they both grow 
approximately according to the scaling ir~ /ir + « (N/Z) 2 
predicted by the A(1232) resonance model but fall far be- 
low the FOPI data. Our calculations with varying values 
of x indicate that a strong symmetry potential is at work 
at this energy as one expects. While not perfectly re- 
produced by our calculations even with x = 1, the FOPI 
data suggest unambiguously that the E sym (p) is rather 
soft at supra-saturation densities compared to the APR 
prediction. 



Since the Au+Au collision is most sensitive to the 
Esy m (p) among the reactions considered here, we now 
turn to the excitation functions of the pion yield and the 
7r~ /tt + ratio for the most central Au+Au reactions. Both 
the pion yield and the tt~ / 7r + ratio do not change notice- 
ably within error bars with the fine sub-division of the 
impact parameter within bo < 0.15. We thus compare 
simply calculations for head-on collisions with the FOPI 
data for the most central reactions. Figure [3] displays 
the excitation function of the pion multiplicity per par- 
ticipant M^/A part . In order to compare with the data 
directly, the total pion multiplicity is obtained from the 
charged pions only by using 1.5 x (M n - + M„+) and the 
number of participants is calculated from 0.9 x A sys where 
A sys is the total mass of the colliding system as done in 
the data analysis [2{|. It is seen that the results of the 
calculations arc in reasonably good agreement with the 
available data. We notice here that, unlike the ir~ /ir + ra- 
tio which has the advantage of reducing significantly not 
only systematic errors but also effects of isoscalar nuclear 
potentials, the pion yield also depends appreciably on the 
EOS of symmetric nuclear matter [29l. [3o|. The inset il- 
lustrates the time evolution of the central density with 
the maximum value varying from about 2.2po to 3.5/Oo 
for the beam energy from 0.25 to 1.2 AGeV. 
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FIG. 4: (Color online) Excitation function of the -k~ /-k + ratio 
in central Au+Au collisions calculated with the IBUU04 in 
comparison with the FOPI data and the IQMD prediction. 
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FIG. 3: (Color online) Excitation function of the pion multi- 
plicity per participant M n /A pa rt in the most central Au+Au 
collisions. The inset shows the time evolution of the central 
density. 



Shown in Fig. [4] are the excitation functions of the 
7r~/7r+ ratio calculated with the IBUU04 and the IQMD 
in comparison with the FOPI data. First of all, the 
decreasing trend of the n~ /n + ratio with the increas- 
ing beam energy is well reproduced by all calculations. 
Most interestingly, IBUU04 calculations with x = 1 
can best describe the FOPI data over the whole energy 
range. Moreover, the it~ /it + ratio is more sensitive to 
the E sym (p) at lower beam energies as expected. 

Putting together all available information wc infer that 
the E sym (p) reaches a maximum somewhere between 
po and 2po before it starts decreasing at higher densi- 
ties. This indicates the importance of the development 
at NSCL/MSU 0] and RIKEN [47| of Time Projection 
Chambers (TPC) to measure ir~ /ir + ratios in heavy- ion 
reactions induced by radioactive beams from about 120 
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to 400 AMcV. These measurements around the E^ h can 
help constrain the E sym (p) at densities from about p to 
2pQ. Moreover, the neutron-proton differential flow mea- 
surements planned at SIS/GSl[48| and possibly also at 
the future FA1R/GS1 will not only test our findings here 
but may also enable the exploration of the E Bym (p) at 
even higher densities. Generally speaking, there are lim- 
itations on how much one can learn reliably about the 
high density E sym (p) from high energy heavy- ion col- 
lisions. For instance, to reach higher baryon densities 
by using more energetic beams, more baryon and meson 
resonances are excited. Moreover, with the increase of 
the beam energy, the collision becomes more transpar- 
ent Therefore, the hot and dense system formed is 
further away from the pure nucleonic matter in thermo- 
dynamical equilibrium. It thus becomes more difficult 
to model and relate properties of this system with the 
E sy m{p) of nucleonic matter. In this regard, neutron 
stars are unique laboratories for probing the E sym (p) of 
cold high density nucleonic matter. On the other hand, 
astrophysical studies of neutron stars have their own lim- 
itations and challenges. It is thus important to do cross 
checking between the two fields. If our conclusion on 
the E sym (p) is confirmed by more experimental and the- 
oretical studies, it not only posts a serious challenge to 



some nuclear many-body theories but also has important 
implications on several critical issues in nuclear astro- 
physics, such as, the cooling of proto-ncutron stars, the 
possible formation of polarons due to the isospin sepa- 
ration instability [ToL [la] . the possible formation of quark 
droplets and hvperons[2(| in the core of neutron stars. 

In summary, circumstantial evidence for a rather soft 
nuclear symmetry energy E sym (p) at p > 2po compared 
to the APR prediction was obtained from analyzing the 
recent FOPI data on the ir~ /tt + ratio using the IBUU04 
transport model. Further experimental and theoretical 
confirmations of our findings are important. 
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